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Abstract 

Initially, we derive a nonlinear integral equation for the vacuum counting func- 
tion of the spin 1/2-XYZ chain in the disordered regime, thus paralleling similar 
results by Kliimper pQ, achieved through a different technique in the antiferro- 
electric regime. In terms of the counting function we obtain the usual physical 
quantities, like the energy and the transfer matrix (eigenvalues). Then, we in- 
troduce a double scaling limit which appears to describe the sine-Gordon theory 
on cylindrical geometry, so generalising famous results in the plane by Luther [5j 
and Johnson et al. [3]. Furthermore, after extending the nonlinear integral equa- 
tion to excitations, we derive scattering amplitudes involving solitons/antisolitons 
first, and bound states later. The latter case comes out as manifestly related to 
the Deformed Virasoro Algebra of Shiraishi et al. [I]. Although this nonlinear 
integral equations framework was contrived to deal with finite geometries, we 
prove it to be effective for discovering or rediscovering S-matrices. As a particu- 
lar example, we prove that this unique model furnishes explicitly two S-matrices, 
proposed respectively by Zamolodchikov [S| and Lukyanov-Mussardo-Penati [HlITj 
as plausible scattering description of unknown integrable field theories. 
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1 Introduction 



By definition, two dimensional (quantum) integrable models possess as many (possibly 
infinite) commuting independent (conserved) quantities as the degrees of freedom. In 
the quantum scenario and especially when the number of degrees of freedom is infi- 
nite, integrability does not still guarantee the access to exact information through a 
simple and standard way. Nevertheless, under propitious circumstances simultaneous 
eigenvectors and corresponding eigenvalues may be obtained by using various methods, 
among which the Algebraic Bethe Ansatz is a very powerful and customary technique. 
One reason of its great success may be tracked in the eigenvector construction, though 
this reveals itself rather cumbersome (but not available in many other Bethe Ansatz 
versions). Another motive may be its huge range of applicability, in particular within 
the spin chain world. As well-known, this may be thought of as equivalent to classical 
statistical lattice models (in two dimensions), which often describe important generali- 
sations of scaling quantum field theories. Moreover, although the vocabulary from spin 
chains to (2D) field theories has still missing items and conjectures in its pages, it seems 
to have recently gained a new section on spin chains hamiltonians as mixing matrices in 
four dimensional (super) Yang-Mills theories ([BUS!, |X0] and the development thereof 
originated) . 

From the physical point of view, one of the most important effects involves how 
the properties of the system - and in particular the eigenvalues of the commuting 
observables - vary with its spatial dimension, i.e. the number of sites in the spin chain 
case. Besides the relevance of finite size effects in statistical mechanics and condensed 
matter (e.g. JJJ and references therein), the scale of coupling (and then of energy) also 
seems to be tuned by the spin chain length in the aforementioned 4D field theories of 
strong interaction. 

As for finite size effects, the Non-Linear Integral Equation (NLIE) description - first 
introduced in ^2] for the conformal vacuum and then derived for an off-critical vacuum 
in by other ways - turned out to be an efficient tool in order to explore scaling 
properties of conserved charges. Since [Hj, still concerning features of the vacuum, and 
[To] , regarding excited states, a number of articles was devoted to the analysis of and 
through a NLIE and mainly follows the route pioneered by Destri and de Vega ^3] (cfr. 
the Hungarian lectures jTHj for an overview). In this way (which will be ours too), the 
NLIE stems directly from the Bethe equations and characterises a quantum state by 
means of a single integral equation in the complex plane. The NLIE has been widely 
studied for integrable models described by trigonometric-type Bethe equations: for 
instance, the 1/2- XXZ spin chain ^7], the inhomogeneous 1/2-XXZ and sine-Gordon 
field theory (ground state in [Tl] . excited states in [T3]) and the quantum (m)KdV-sG 
theory [JSj . Here instead we wish to understand better a less studied, but more general 
set-up: the elliptic one. In particular, we choose the 1/2-XYZ spin chain as prototype of 
our investigation, since it is the direct generalization of the basic trigonometric models 
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and yet not so much complicated to prevent a detailed and profundus analysis. Its 
hamiltonian may also have serious chances to represent a mixing matrix in some gauge 
theory. 

In addition, we have at least another motivation to study elliptic theories and es- 
pecially elliptic spin chains. Again this comes out partially from 2D integrable field 
theories. The latter are often studied by starting from a scattering S-matrix, which 
replaces somehow the Lagrangian as definition of the theory on the plane (the most sig- 
nificant example being the sine-Gordon study pioneered by both Zamolodchikov [TTij). 
Very often both S-matrix and Lagrangian (or the specific perturbation of a conformal 
field theory) are well identified and tied together. Therefore, the S-matrix has naturally 
become the definition of a field theory, provided it verifies all the field theory axioms 
and integrability. Nevertheless, the field theory counterpart of a specific S-matrix is 
sometimes not clear in the Lagrangian language and the assumption of the S-matrix as 
starting point should be only a useful working hypothesis towards further investigations 
and identifications. The elliptic case is indeed a good example of this phenomenon with 
its beginning [S] coeval with the trigonometric relative [T^j, which is crystal clear and 
paradigmatic since then. On the contrary, the whole elliptic scenario is still contro- 
versial and very subtle, although some proposals were recently supported [2D]- m this 
article we want to modify the current perspective on the correspondence problem and 
prove that the 8- vertex Hamiltonian is responsible for all the known elliptic S-matrices 
and gives naturally a unitary explanation of their appearance: this effort is willing to 
give a global view on a fragmented matter. Of course, we might also interpret this as 
a step towards a field theory description or as an interesting tool to generate elliptic 
S-matrices (e.g. by raising up the chain spin). In any case, scattering amplitudes may 
be naturally obtained from the finite size set-up, even though these are clearly defined 
when the volume is infinite. Moreover, as first pointed out in their derivation 
requires some information about the excited states. More specifically, if the connection 
with the kernel of the NLIE may appear a promising gloss in ^Hj and a comprehensive 
elaboration in [21] as regards the sine-Gordon field theory, it can be regarded here as a 
very profitable, predictive and general tool of investigation and analysis. 

Concerning the content and organisation of this article, we address the problem of 
writing a NLIE satisfied by the counting function of the spin 1/2-XYZ chain in the dis- 
ordered regime, by following the more versatile route of Destri and de Vega [T3*| IT^ ITK] 
(Section 2 about the vacuum and beginning of Section 6 about excited states). As for 
the vacuum state, Kliimper wrote down a similar NLIE for the spin 1/2-XYZ chain in 
the antiferroelectric regime and by means of the technique initiated in j!2j . In Section 

3 exact expressions of the eigenvalues of the transfer matrix (and in particular of the 
energy) are given as nonlinear functionals of the counting function. Both consist of two 
terms: a contribution proportional to the size and a finite size correction to it. Section 

4 is devoted to the trigonometric limit towards the conformal (massless) XXZ chain, 
both on the NLIE and on the transfer matrix eigenvalues. Separately (Section 5), the 
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(trigonometric, but massive) double scaling limit - which yields the sine- Gordon field 
theory on a cylinder - is performed and then compared to the results of [TH] • In Section 
6, by considering the first excitations upon the vacuum state, we are able to derive 
the scattering S-matrix of the soliton/antisoliton sector (in the repulsive regime): this 
matrix turns out to have an elliptic form and coincides with that proposed by general 
principles (field theory axioms; integrability: factorisation and Yang-Baxter relation) by 
A. B. Zamolodchikov some time ago |SJ. Furthermore, we compute the scattering factor 
of the lightest soliton-antisoliton bound state in the attractive regime: on the contrary, 
this manifestly coincides with the structure function of the Deformed Virasoro Algebra 
(DVA) by Shiraishi-Kubo-Awata-Odake 0] and hence reformulates in XYZ variables 
the factor hinted by Lukyanov [6 . Then, we thought of this factor as that describing 
the scattering of the fundamental elliptic scalar particle, and therefore as a suitable 
candidate to describe an elliptic deformation of the sinh-Gordon theory And indeed, 
after another mapping of variables, it coincides with the starting definition adopted 
more recently by Mussardo and Penati [7j, although the cumbersome necessary algebra 
has made this coincidence likely unnoted (cfr. also [20J). Therefore, we find an unified 
arrangement for both previously known elliptic S-matrices, though the underlying the- 
ory is not properly a field theory: on one side the elliptic Zamolodchikov S-matrix, on 
the other the DVA or Lukyanov-Mussardo-Penati factor. Eventually, some conclusions 
and many perspectives come to mind and part of both is outlined in Section 7. 



2 The lattice theory 

The spin 1/2-XYZ model with periodic boundary conditions is a (lattice) spin chain 
with hamiltonian written in terms of Pauli matrices a x,y,z 

1 - 

n = -- X)( J*« + i + Jy«+i + Jz«+i) • (2-1) 

n=l 

Here N is the number of lattice sites and because of the periodicity the site N + 1 
is identified with the site 1. The three (real) coupling constants J x , J y and J z may 
be reparametrised (up to an overall constant) in terms of elliptic functions (as for the 
notations on elliptic functions, we refer to [22.). In fact, after introducing the (complex) 
elliptic nome q (\q\ < 1), we may define the modulus k (and the complementary modulus 
k' = Vl — k 2 ) and the associated complete elliptic integral of the first kind K: 
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We also introduce the parameter K' such that the elliptic nome reads 

q = . (2.3) 

The theta-functions of nome q, 6 a b (u; i^-), a, b = 0, 1, may now be defined as 

dooUi^J = (-g e 2 -;g 2 )(-ge~ 2 -;g 2 )(g 2 ;g 2 )=^(v™ ; z^ 

Ooi U i^) = (qe 2l7TU ; g 2 )(ge" 2 "™; g 2 )(g 2 ; q 2 ) = 9, L U] i^) , (2.4) 

010 ( u ^) = 2gicos7r M (-g 2 e 2 ^;g 2 )(-g 2 e- 2 ^;g 2 )(g 2 ;g 2 ) = 6 2 (ttu;^ 

011 = - 2 ^sin7r M (g 2 e 2i ™;g 2 )(g 2 e- 2i ™;g 2 )(g 2 ;g 2 ) = -B x (ttu;^ 
where we have introduced a shorthand notation for the infinite products, 

oo 

(x;a) = JJ(1 -xa s ) , (2.5) 

s=0 

and at the furthest right the expressions in terms of the usual Jacobi theta-functions 
9i, i = 1,2,3,4. With the exception of Section 6, we will find more convenient to use 
the alternative functions: 



e ( w ;^J = ^oi f 7^;^ J, ©i J =0oo f 7^;^ ) ■ (2-7) 



k; V2k'"k; ' ~ x \~'~kj " uu ("2K'"k 

The Jacobian elliptic functions are 

" Vke(«;if)' cn rKj Vke(«;f) 



dn «;<— ) = y/g-^Lgl. (2.9) 



In order to simplify notations, from now on we will omit the dependence on the elliptic 
nome, when the elliptic functions have elliptic nome q ()2.3I) . 

In terms of Jacobian elliptic functions the coupling constants J x , J y and J z in (j2.1j) 
are parametrised as 

J x = 1 + k sn 2 2r? , J y = 1 - k sn 2 2r? , J z = cn2r? dn2r/ , (2.10) 

where 77 needs to be real. The eigenvalues of the spin 1/2- XYZ transfer matrix were 
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firstly found by Baxter [231 by using his auxiliary matrix (Q-operator) technique. Af- 
terwards, in the paper [21] Takhtadjan and Faddeev wrote down also the eigenvectors 
(at least formally), making use of the Algebraic Bethe Ansatz. In the present paper, 
we want to restrict our analysis within the disordered regime, i.e. 

< q < 1 , 

0<r]<K, (2.11) 

which in particular entail K > and K' > 0. With this choice of parameters, we differ 
from pQ, where a NLIE was written (by means of a method far from ours, though) in 
the antiferroelectric regime: < q < 1 and < if] < K'. The Bethe equations obtained 
in [231 EH are: 



H{iotj + r])Q(iaj + rj) 
H(iacj — T])Q(iaj — rj) 



N n 



-i-r H(ictj - ia k + 2r])Q(iaj - ia k + 2rj) 
H(iaj — ic*k — 2rj)Q(iotj — iotp, — 2rj) 



(2.12) 

for all j = 1, ...,n, where v is an integer. Equations ()2.12)) are valid when 

m 1 r] = 2m 2 K, (2.13) 

where m x and m 2 are integers, and when 2n = N (mod rrii). The corresponding 
eigenvalues of the transfer matrix are 

A N (a) = e^e(0) N H( ia + r,)»G(ia + y) N f] H ^~ ia ± 2 ^fi" Z * a + 2 ^ 

J-X hi in/ . — i rv 1 M f i rv ■ — ■? rv\ 



] tl (totj — ia)& — la) 

(2.14) 

H{ia — i(Xj + 2r])Q(ia — iaj + 2rj) 



if (za — ?o;y)B(zQ; — %a 
In writing (|2.121 12~TH) we used notations by |24j . 

2.1 Nonlinear integral equation for the vacuum 

We now want to study the Bethe state with lowest energy, i.e. the vacuum. From 
[2U 123 we know that such a state is given by all real roots (ctj G M.) of the Bethe 
equations enjoying these additional properties 

N K' K' . 

u = 0, n = —, -—<«,<—. (2.15) 

Using symmetry properties of elliptic functions, we rewrite equations (J2.12|) as 



H{ia.j + r))Q(i(Xj + rj) 
H(—ictj + r))Q(—iaj + rj) 



N N/2 



TT H(iaj — iak + 2r/)Q(iaj — iak + 2rf) 
H(—iaj + iak + 2r])Q(—iaj + iak + 2rf) 

(2.16) 
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We then define the function 

^ , H(£-ix)G(£-ix) _ m , 

^•fl = ' h g (I + fa )e(| + i,r {eR - (2 ' 17) 

which has branch points occurring at the points x = x r>s of the complex plane such that 
|Rex rjS | = rK' , \lmx ryS \ = £ - 2sK , r, s G Z. (2.18) 
Therefore <f>(x,£) is analytic for x in a strip around the real axis defined by the condition 

\lmx\ <min{£, |£-2K|}. (2.19) 
Having introduced the function ()2.17j) . we may write the Bethe equations as 

N/2 

iN<f>{a h ri) = ln(-l) 1+ ^" + i ^ <f>(aj - a k , 2r]) . (2.20) 

k=i 

We define the counting function 

N/2 

Z N {x) = N<f>(x, v)~J2 ^ x ~ a ^ 2 ^ ' ( 2 ' 21 ) 
fc=i 

which is analytic, as a consequence of f)2.19j) . in the region (containing the real axis) 

2 

\lmx\ < 7] , if < rj < -K , 

3 

2 

\lmx\ < 2K - 2r] , if -K<r^<K. (2.22) 

3 

In terms of the counting function, the Bethe roots are identified by the condition 

Z N ( aj ) = n (21, + 1 + jj , I 3 eZ. (2.23) 
For the sake of simplicity, we restrict our analysis to the case N e 4N, so that 

Now, fundamental property of the vacuum roots, ctj, is not to allow any missing root 
(hole) in equation ()2.23|) and therefore entail a simple sum property about a function 
g{x) analytic around the real axis 

N/2 k!_ k!_ 



2iriJ2g{a k ) = - I ' dxg'(x-ie) In [l + e iZN{x - ie) ] - I dxg'{x+it) In [l + e iZN{x+ie) ] 

k=l " 2 ^\ 

(2.25) 



In (I2.25J) e > is a parameter small enough to keep the integration within the analyticity 
domain 1 (prime means, as usual, derivation). In the limit e — > this equation can be 
rearranged as 



N/2 Kf_ K 



2nJ2g( a k) = -J * ; dxg'(x)Z N (x) + 2 J ^ dxg'(x)lm\n [l + c iZn{x+i0) ] . (2.26) 



fc=l 2 2 



Thanks to the analyticity property (|2.22|) . we may therefore rewrite the vacuum count- 
ing function (|2.21j) in a useful form 



K K' 



Z N (x) = N(f>(x,T])- I * %'(x-y,2 V )Z N (y)+ [ 2 %'(x-y, 2 V )Im In [l + e «»(v+«»] 



2 2 



(2.27) 

where 4>'(x, tj) is the x-derivative. An inspection of the transformation properties of the 
elliptic functions 



H(u + iK!) = iq ™k8(u) 
Q(u + zK') = iq~^e~^H(u) 



shows that <f>(x,£) and, consequently, Z N (x) are quasiperiodic in their regions of ana- 
lyticity (j2.19j) and (|2.22jl respectively, with a real quasiperiod K': 

4>(x + K', £) - (f)(x, = ^U-jQ (x in strip (jni) , (2.28) 
Zjv(a; + K') - Z N (x) = nN (x in strip . (2.29) 

Obviously, the derivatives of <f>'(x,£) and ^(x) are periodic (with period K'). A good 
way to rearrange (|2.27|) is to introduce the Fourier coefficient 

r— 

f(n) = / 2 dxf(x)e 2l W , (2.30) 

for a (quasi)periodic function /(#). In terms of the coefficients f(n) the (quasi)periodic 
function f(x) is expressed as 



-, +oo 

W = ^ E /("K*^ , (2-31) 



K 

n=— oo 

within the principal interval —^-<x<¥f. If we consider a periodic function / and 
(in general) a quasiperiodic function g, it is well known that their periodic convolution, 

(f®g){x)= [ 2 dyf(x-y)g(y), (2.32) 



1 We are implicitly assuming the reasonable hypothesis that e may also be made small enough not 
to allow spurious solutions of 12.24JI within the integration contour. 
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has Fourier coefficients given by the product 

(7®ff)(n) = f(n)g(n) . (2.33) 

Bearing this property in mind, we introduce the shorter notation 

L N (x) = Imln [1 + e lZN{x+l0) ] , (2.34) 

and may easily prove that the Fourier coefficients of all terms in relation ()2.27|) satisfy 
the relation 

Z N {n) = N4>{n, r,) - ^-$'{n, 2r,)Z N {n) + -&(n, 2 V )L N (n) . (2.35) 

271 IT 

And this immediately entails 

Z N (n) = N — — + 2 zn x -L N (n) . (2.36) 

l + ^(n,277) l + JL^( n ,277) 

Being the functions Z^(x), <fi(x) and L^{x) all odd, their Fourier coefficients are van- 
ishing for n = (equation ()2.36|) is trivially satisfied), and hence we may drop out the 
zero mode in the series expansion. Therefore for the Fourier series we obtain 

Z N (x) = NF(x) + 2 / 2 dyG(x - y)Imln [l + e iZN{y+i0) ] , (2.37) 
where we have defined the forcing term 



Fix) 



— Y e" 2 ^, (2.38) 

K' 1 + i^(n, 2rj) 



and the operator kernel 



+ 0O 

G(x) = - V - **YZ~" > (2-39) 



K, ^l + i^(".2i ? ) r 



where we have decided to remove the zero mode. 

Equation (|2.37|) is the Non-Linear Integral Equation describing the vacuum of the 
spin 1/2-XYZ chain in the disordered regime (J2.11|) . Even though it has been derived 
supposing x real, it is valid in the strip (J2.22j) thanks to analytic continuation, which 
would fail at some points at the border of the strip (with real parts given by the first of 
(|2.18|l ). Now, we need to compute explicitly the Fourier series for the functions involved 
in dUHZD- 
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2.2 Calculation of Fourier coefficients 

First, we recall the Fourier coefficient of 4>'(x, 2rj) from the results of Appendix A. Since 
x is in region (j2.22|) . we use the first of (|A.14|) or the second of (|A.15|) to obtain in 
either case 

sinh MK-^)* 

nnM=2* -pi, - (2.40) 
Now, we are ready to compute the coefficients 

0(n, 7])= I dx(p(x,T])e 2in ^ L . (2.41) 

The zero mode 0(0, 77) = is simply given by (p(x, rj) = — 0(— x, 77), while the others are 
given upon integrating by parts and exploiting the quasi-periodicity ()2.28|) : 



4>'(n,rj) = / dx(j)'(x,r))e 



1) 



K' \ ( YJ 



2 , , , N 2imr 

-e ~K r 



, dx(j)(x,r]) — 



= 271 V ~ k) cos nn ~ ~K^ (n ' *>) ■ (2 - 42) 

Plugging (J2.40)) into, we are given the required expression 

0(0,77) = 0, 

iK' f sinh 2ra(K ~ ??)7r / n \ ) 
d>(n,r,) = — I ~p|b -(l-^jcosnTrj, n^O, (2.43) 

which entails, with quasi-periodicity at hand, 

+00 2i n7T ? • 1 2n(K— n)ir 

sr^ e ^ sinh — v y , - / 77 \ irx . . 

0(x,77)=iV ^ — + 2M--M— . 2.44 

y " ^ n sinh^ \ KJ K> K ' 

n= — 00 xv 

And we may re-write this series in a more compact form as 

+2; sin Us* sinh 2 "(* 7 » 
«'■■>)= E / "pfi ■ (2-45) 

n=— 00 xv 

whose convergency domain is (like in (|2.44j) ) |Imx| < 77, containing the strip (|2.22j) . 
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2.3 Forcing and kernel functions 

Given the preceding expressions of 0(n, 77) ()2.43j) and of 0'(n, 77) (|2.40)l . we may explicitly 
mould the forcing term (|2.H8jl into 

/ x i sinh 2 "^7'^ 7r — (l — IL) cos nir sinh 

F(x) = V ^ i o ^ ^ — e~ "k 7 " . 2.46 

K) ^ n sinh + sinh 2^-2^ V ) 

n= — 00 J\. Jv 

From ()2.39|) we read off the Fourier coefficient 

a (B)= ^"' 2 "> , ( 2 .47) 
1 + ^(n, 2,) 

and then simplify this by means of (|2.40p . 



sinh 



2n(K-2-q)iT 



G(n) = , K ' , (2.48) 

2 sinh cosh ^ V ' 



whence to obtain the final expression for the kernel 

G{x) = ±y — e~ 2 ^ . (2.49) 

K J K' ^ 2 sinh 2 "^7^ cosh^ 1 ' 

71— — CO J£' K 

Remark. In the case 77 = K/2 the kernel function G(x) vanishes. Therefore the 
solution of the NLIE is trivial: 

+ °° sin 2n7TX f K \ / K \ 

Z N (x) = N > = iVam 2x; z— = iVarcsinsn 2x; 7— , (2.50) 

v ' ^ 2ncoshi# V KV V K/ 



K 



where the elliptic functions have nome q' = e -7 ^. This is not a surprise as this case 
reduces to the free-fermion point 77 = 7r/4 in the trigonometric (q — > 0, cfr. Section 4) 
limit. 



3 Eigenvalues of the transfer matrix 

It is the main aim of this Section to write down the eigenvalue of the transfer matrix 
(12.14)) on the vacuum state in terms of the solution of the (vacuum) Non-Linear Integral 
Equation ([2.37)1 . Although this might seem a priori a complication, it has at least two 
main advantages. First, instead of solving a big number of transcendental (Bethe) equa- 
tions, one should solve a NLIE: this makes numerical computations and approximations 
much simpler. Second, in the following expressions the bulk terms (proportional to the 
size N of the system) are clearly separated from their finite size corrections. And as it 
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will be clear later on, the finite size corrections and properties can be singled out in the 
limit N — > oo. 

Again, we start from relation (j2.26|) which expresses the sum on the vacuum Bethe 
roots of an arbitrary analytic function g. Using the NLIE, ()2.26j) may be written as 
follows 

N/2 K' 
Z-y N f~ 

2_^g{uk) = ~y I K , dxg{x)F{x) + 

k=l ° 2 

+ / 2 —g'{x) I 2 dy[5{x~y)-G{x-y)]lmln\l + e iZNiy+l0) ] . (3.1) 

/ K' 7T / K' 

J a" J 2" 

In terms of the Fourier coefficients this relation reads as well 

N/2 +oo 

E«N = -^; £ Hn)g'{-n) + 

k=l n=—oo 

+ r I!?.!!, y ^(n)[l-G(n)]e-Tllml n ri + e ^(^)l . (3.2) 

J _KJ_ 71 K' 



n=—oo 



Therefore, once the solution of ()2.37|) is found, the relation ()3.2|) may be used in order 
to compute (very often numerically) the eigenvalues of observables on the vacuum. We 
remark that the first term in ()3.2)1 is apparently proportional to N . Therefore it gives 
usually the leading order in the N oo limit. In fact, the second term in (|3.2j) usually 
provides the sub leading corrections in that limit. It is convenient to define the function 
constructed upon g(x) 

J ^) = k7 E 9'{n)[l - G(n))e- 2J ^ , (3.3) 

n=— oo 

which appears explicitly in ()3.2|) . Here we decide to add to the definition of G(x) its 
zero mode, since its inclusion does not change the value of the integral in ()3.2|) and 
makes at the same time formulae more compact. In particular, we want to apply the 
result ()3.2j) to the vacuum eigenvalues of the transfer matrix ()2.14j) . In this respect, we 
may introduce the decomposition 

A*(«) = BiOfHiia + „)*6(ia + v) N TT ^V,™ + ^nfi" " ~ + 

j\ M \ %a j ~ ioi.)<d\ia.j — ia) 

+ B(0)» H(ia - V ) N e(ia - V f ff ^ ~ ^ + 2 ^' a - ^ + ^ = 

H(ia — iaj)Q(ia — iaij) 

= A+(a)+A^(a). (3.4) 
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Let us first concentrate on lnA^(a). We have that 

N/2 

lnA+(a) = iVlne(O) + NlnH(ia + rj) + Nlne{ia + rj) + ^7 + (a fe ,a) , (3.5) 

k=l 

where 

. , , H(ix — ia + 2r])&(ix — ia + 2rj) 

-y+{x,a ) = ln — — — — . 3.6 

H [ix — ia)<o{ix — la) 

Because of the periodicity property 

lnA+(a + 2iK) = lnA+(a), (3.7) 
we can restrict the complex a to the strip 

- 2r] < Ima <2K-2rj. (3.8) 
Upon comparing (|3.6|) with (|2.17|) . we see that 

7+(x, a) = —i(j){x + zK — in — a, K — 77) . (3.9) 
Recalling (jA.13|) we easily obtain the Fourier coefficients 

— 

7',m, a) = / ax— r y+(x, a)e , (3.10) 
J-^ dx 

in a form depending on the imaginary part of a 



y J n , a) = 2m ^ — e ^(™-«) if - 2r/ < Ima < , 

' +\ > J sinh 2nK7r ' 



K' 



YJn,a) = -2m -^ e -^( K +^-") , if < Ima < 2K - 2n 



sinh ^7 



(3.11; 



In this case, the function J 7+ (|3.3|) equals 

+ OO 2injr r . ■ N 

. . ITT ^-^ e ~K>-( a + lr l X > 

Jj + (x,a) = — 2^ CQoh 2n^ > if -2r? < Ima < 0, 

n=— oo K' 

(3.12) 

J 7 (x, a) = -g- V 9 5 e^ 8t * HK -'), ifO<Ima<2K-2r / . 

K' n f^ sinh Hlifczk cosh iHZZL 

As a consequence, from relation ()2.26|) and from the expression for F(n) ()2.46|) we are 
finally given the two cumbersome expressions: 
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' if —2r) < Ima < 0: 

= N In 6(0) + 



+ 



lnA+(a) = N 



-\-o(j 
n=— oo 



n 



N In H(ia + rj) + N In ®(ia + rj) + 
tt, (1 — |£jcoswr 



sinh g|ggE (i _ |) cos 

2 cosh ^ sinh ^ 2 cosh § 



+ 



+ 



2 z ^-^ e k' ^ a+iV x > r 

cfe — > ^ Imln 1 + e 

k, K ' ^ cosh^ L 

2 n=— oo K' 



(3.13) 



\z N - 



•(aH-tO)j 



if < Ima < 2K - 2r?: 

lnA+(a) = JVlne(O) + N\nH(ia + r)) + 



n= — oo 



n 



- N\nQ(ia + rj) - 
sinh 2MZL (i _ i) CO s nvr s inh ^ 

2 cosh ^ sinh ^ ~ 2 cosh ^ sinh 2w( *T^ 

k' K' ^ sinh 2 <k-^ cosh ^ L J 

r n=-oo &llm K' ousu K' 

For what concerns lnA^(a), we have that 



(3.14) 



In A N (a) 
where 



N/2 

= N In 9(0) + NlnH(ia - rj) + N In B(ia - rj) + ^7_(a fc , a) , (3.15) 

k=l 



, , , — zx + 2n)Q(ia — ix + 2?i) 

7_(x, a) = In — — — — . 

H [la — ix)<a{ia — ix 



(3.16) 
for 



We notice that j-(x, a) = —■y + (x,a — 2%r\). By using this link, the expression 
lnA^(a) may be easily obtained from (|3.13l I3~T4^) . 

3.1 The energy 

It may be of some interest to write down the vacuum eigenvalue of the Hamiltonian 
(|2.1|) in a form which separates the term proportional to N to its finite size correction. 
From [21] we learn that in general we may extract the eigenvalues of the Hamiltonian 
(of the spin 1/2- XYZ chain) from the eigenvalues of the transfer matrix as 



En = isn2r] —— In An (a) 
da 

a=—ir] 

which in the notations of last subsection takes on the form 

d 



En = isn2i] — In A N (a) 
da 



(3.17) 



(3.18) 



a=— irj 
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In the particular case of the vacuum, result ()3.13|) yields the wanted expression 



E 



(vac) 
N 



-Nsn2r] 



H'{2rj) &{2rj) 
_H{2rj) + &(2ri) 



+ 



+ isn2 V \ -N ^ — 



sinh 



2n(K-r])w 
K 7 



cosn7r 



Ki K' 

2 n=— oo 



2 cosh ^ sinh ^ 2 cosh ^ 

]}■ 



+ (3.19) 



i - , l \ -v 2ot7t e k r 7 , 



As easily follows from this explicit expression, the last term gives the finite size correc- 
tions to the first two terms in the limit N — ► oo. 



4 Trigonometric limit (i.e. XXZ chain) 

In order to have a check on the validity of our results and to prepare the ground for the 
more refined limit of next Section, it is important to explore the trigonometric limit in 
which the spin 1/2-XYZ chain reduces to the spin 1/2- XXZ chain (q — > or J x = J y 
in (|2.1|l ). We expect to reproduce the results of papers [T71 [T3] • As the trigonometric 
limit is also expressed by the limit K' — > +oo (therefore K — > 7r/2), a Fourier sum can 
be replaced by an integral, according to the prescription 



1 +oc / n \ r +oc 



In this Section we aim at writing down the trigonometric limit of the vacuum NLIE 
()2.37)l and of the vacuum eigenvalues of the transfer matrix ()3.13| 13.14)1 and of the energy 
(J3.19)) . With reference to (|2.11|1 . the range of variation of rj becomes < rj < tt/2. By 
applying prescription (|4.1j) to ()2.49|) . we obtain an integral expression for the kernel 
function G(x) 



/+oo 
dpe 
-oo 



-2ipTTX 



sinh [2p (f - 2rj) tt] 



x 2 sinh \2p f~ — rA it] cosh 2pr]ii 



" e i P x sinhp(f-2^) ^ 
s 2tt 2 sinhp f~ — rj) cosh.prj 

As a check we may notice that the inclusion of the zero mode of G is irrelevant because 

±6(0) = -L*-H (4.3) 
K' 2K' K — T] K ' 

is vanishing in that limit. Analogously the trigonometric limit of 

. +oo 

<t>'{xM = ^,Y. e- 2i ^$'(n,2 V ) (4.4) 



15 



may be easily computed starting from ()2.40|) : 

[ + °° . sinhp (f - 2ri) 

^^L ipe (4 ' 5) 

For what concerns the forcing term F(x), in the trigonometric limit the oscillating term 
in ()2.46|) gives no contribution and leaves us with 

1 f +ca dp sinpx 1 / mx tt\ nx 

Fix) — > — I = - In tan 1 — = arctansmh — . (4-6) 

2 J_ OQ p cosh pr] i \ Ar) 4 / 2r\ 

Therefore, in the trigonometric limit the NLIE ()2 . 3 7|) takes on the form 

Z N (x) = Narctansinh — + 2 / dyG sG {x - y)Imln [l + e lZN{y+l0) ] , (4.7) 



which is the Non-Linear Integral Equation for the spin 1/2-XXZ model in the massless 
antiferromagnetic regime (i.e. J x — J y — 1, J z — cos2rj). Moreover, we may straight- 
forwardly perform the limit on the vacuum eigenvalue of the transfer matrix, or better 
on (j3.13| l3~l"4j) . with the outcome 

• if — 2r] < Ima < 0: 



: +QO dp ■ , siahp(^-ri) 
lnA+ a = N / ± e ~^ P{ \ 'I + 
' p 2 coshp?] smh ^ 



(4- 



/+oo 1 
dx -. -1mm [l + e iZN{x+i0) ] 
2rj smh^(a-x) 

if < Ima < 7i — 2rj: 

In A+(a) = -N cte-wa+PV-pZ sinh ^ 



x p 2 cosh pr] sinh 2 



^ - / .. ' ■ o, 1, ;„i, VJL 

iz N (x+ioy 



(4.9) 



- f °° dx f °° ± ; Sinh ^ e-(^f-)l m ln [1 + e 

J-oo J-oo 27r sinh p (| - rj) cosh prj 

Eventually, the energy (j3.19j) takes on the limit expression 

E -> -N cos 277 + sin 2r] / dp— V2 



sinhp7r 2 cosh 2pr]7T 



(4.10) 



/"+oo sinh 

- — sin 2ri / dx ^-Im In [l + e iZN{x+l0) ] 

4r/ 2 '7-eo cosh 2 ^ L J 



2,, 

And indeed these formulae coincide with those given in f7|EI] as for the spin 1/2-XXZ 
model. 



16 



5 A double scaling limit: cylinder sine-Gordon 



In this Section we want to rewrite the NLIE ()2.37j) in a particular limit such that it 
will describe the sine-Gordon field theory on a cylinder. This procedure extends the 
celebrated results of [2] and [3] which concern how a peculiar infinite length 2 limit of 
the XYZ reproduces sine-Gordon on the plane. We first notice that expression ()2.46|) 
for the forcing term may be separated as 

"°° +^ i (l-A)sinh 2 ^ 



NF(x) =NT — 6 N T — ll "^ Mnnir e -«gW ( 5 1) 

{ ) ^ 2ncosh^ ^ 2n sirih M^ cosh ^L e " ^ 

n— — oo xv n— — oc Jv K 



Then, we shift the variables x and y in (|2.37|) as x = x' + 6 and y = y f + Q, assuming 
6 x-, ^/-independent and positive: 



Z N {x' + 0) = iVF(x' + 0) + 2 / dy'G(x' - y') Imm 



4--e 

i^iv(2/'+e+«0) 



1 + e l 



K' K' 

-— -0<x'< — -6. (5.2) 
2 2 K J 

Now, the double scaling limit is realised by allowing K' — > +oo when N — > oo according 
to 

K' = din DiV , d>0,£>>0, (5.3) 
provided that also diverges as 

= clnCiV, c>0,c<^,C>0. (5.4) 

The reason why this is a double limit (namely the lattice spacing A is going to zero 
too, whereas the length R = NA remains finite) is encapsulated in the coefficients C 
and D and will be clear below. Indeed, since the shift (|5.2|) we had in mind to define a 
modified counting function too: 

Z(x') = lim Z N (x' + 0) . (5.5) 

N^oo 

Of course the range of variation of the new independent variables x' and y' gets more 
involved 

-~ln£>iV -clnCN <x',y' < -InDN - clnCN . (5.6) 

Now, let us derive the limiting value of the forcing term NF(x' + 0): the sums in (J5.1)) 
may be replaced by integrals according to the rule (|4.1j) . i.e. 

iN f +ao dp e - 2i P 7T ( x '+ clnCN ) 



iVF(a/ + 0) = ^ J 
iN f +co dp (l - f ) sinh jot 2 e~ 2il 



p cosh 2prj7T 



2ipn(x'+c In CJV-f In DN) 



2 7-oo P sinh \2p-n f| — 77)] cosh 2^7771 



(5.7) 



2 The length is defined as R = NA, with the lattice spacing A. 
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Because of properties ()5.3J) . (j5.4[) . the first (second) integral of this expression can be 
calculated by closing the integration contour in the lower (upper) p-comp\ex half plane 
and avoiding with a semicircle the singularity at p — 0. In both cases, the leading 
contribution is given by the pole with the smallest modulus. When tt/6 < i] < tt/2, 
these poles are at p = ±i/(4r]) and therefore the total contribution reads 

-1 C7T C7T x' 7T f 2 77 \ 1 I 7T Z' „ (t\ 7T^ 7TC Tvd x' TV 

Nn - 2N 1 ~^C~^e~^ +21 '- N ^ 2 > tan — C^D~^e^ . 5.8 

\ TV J Arj v ' 

The first term comes from the semicircle around p = 0. For technical reasons we need 

2 

to restrict further the domain of r\ to within | < r\ < j, so that tan ^ > 0, and choose 
2r7 , 8ri „ 4 „ 16 / 2ri\ vr 2 . . 

7T 7T mit m Z K Z \ 7T J %] 

where m is a positive constant with the dimension of a mass and R is the lattice length. 
This choice entails 

NF(x' + 0) = Ntx + mR sinh — + o(N°) , (5. 10) 

2r/ 

where o(z) means "order less than z" . Of course, the prescription (|4.1j) produces on 
the kernel function G(x' — y') ()2.49|) the result (|4.2|) as in the trigonometric XXZ limit. 
Therefore, we obtain that the limiting counting function Z(x') ()5.5|) has to satisfy the 
equation 

tit' f +OQ r i 

Z(x') = mRsmh + 2/ dy'G sG (x' - y')lm\n 1 + e lZ(y ' +i0) , (5.11) 

where, since iV G 4N, the constant Ntt has been reabsorbed in a redefinition of ^(x') 
and eventually the interval (|5.6[) has become infinite 

— oo < x',y' < +oo . (5-12) 

With the identification 

' = \ ( x - 9 ■ < 5 - i3) 

the Non-Linear Integral Equation ()5.11|) describes the vacuum of the sine-Gordon field 
theory with coupling constant b 2 and renormalised mass parameter m (i.e. Lagrangian 

2 

£ = \{d(j)) 2 + cos 60) on a (space-time) cylinder with spatial circumference R. 

Remark 1 The Non-Linear Integral Equation (|5.1H1 has been obtained when 7r/6 < 
7] < 7r/4. However, it can be considered without problems in the whole region < 
f] < 7r/2, defining everywhere by analytical continuation the state which minimises the 
energy. 

Remark 2 We remark that in formula (j5.1|) one could have made the choices 
e -2i^{x+±K') Qr | e -2i^(a;+^K') + | e -2if5-(x-±K') in the exponent i a i Q f ^ e j ast term, in 
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order to express the factor cosn7r in (|2.46j) . However, the first choice together with the 
choice x = x' — 9, y = y' — 6, leads to the equality 



7T0C 

NF(x - 6) = —Nn + mi? sinh — + o(N°) 

2,7] 

Instead, the second choice together with the shifts x = x' ± ©, y = y' ±9, gives 

1 ITT 

NF(x ± 6) = ±-Nn + mR sinh + o(N°) 

2 2r\ 

Since the constants —Nn, ±|iW are inessential, we eventually obtain again equation 

Remark 3 With the choices f|5.3l 15.41 15. 9j) . in the limit iV — > oo the parameters J 2 , 
J y , J z (|2.1(J|) behave as follows 

J 2 



(5.14) 



(5.15) 



— — > cos 2r] . 

<J T. 



j X jr.. 



( MR\ 
\4N 



8,, 



J x 8 sin 2 2?] 

where we have introduced a rescaled mass parameter 

M = 



rn 



(l"?)tang 



(5.16) 



(5.17) 



Clearly, the famous scaling limit to the sine-Gordon field theory on the full spatial line 
j2 El is now gained from here by sending R —>■ oo. 

Remark 4 The limit discussed in this section can be applied to the eigenvalues 
()3.131 13.14|) of the transfer matrix. The calculations of the limits of the Z-independent 
terms are less straightforward, but they are carried out following a procedure analogous 
to that used on the forcing term (remembering that due to the redefinition of Z the 
forcing term F(x) has to be replaced by F(x) — it). Here we just give the result: 

• if — 2r] < Ima < 0: 

r 2 



7T 7TCX 

lnA + (a) = — m-Rcotan— cosh — 

4?7 2?7 



(5.18) 



+ 



if < Ima < 7r — 2rj: 



dx 



1 



2r] sinh ^(a — x 



-Imln [l + e ^ + *°)] 



1 



lnA + (a) = —mR j cosh 

sin 2- 

+ °° dp isinhp?] 



71 


f.TT \" 




% a) 


_2r] 


V 2 /_ 



dx 



27rsinhp (f — ^) cosh pi] 



^tpla+irj-i 



(5.19) 

5-*) Imln [i + e ^+ J °)] , 
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which hold when 7r/4 < rj < tt/2. Thanks to a sign change of a, these unveil the 
coincidence with the vacuum eigenvalue of the transfer matrix in the sine- Gordon theory, 
A sG (a), as computed in the second of (HI 

lnA+(a) = lnA+ G (-a) . (5.20) 



6 Scattering theory 

As the spin chain becomes infinitely long, A" — > oo and fixed A, the infra-red (IR) 
Hamiltonian eigenstates become on mass-shell states. Besides the statistical mechanics 
interpretation as thermodynamic limit configurations, these may also be thought of 
as describing the asymptotic particles in the quantum field theory perspective. And 
in the same view the excited state version of the nonlinear (elliptic) integral equation 
()2.37j) . being merely a quantisation rule, gives exactly the scattering amplitudes of 
the asymptotic states provided that the convolution integral is neglected (when A" — > 
oo). Moreover, integrability reveals itself useful in the factorisation of the general 
scattering S-matrix into two particle ones, which eventually encode all the on mass- 
shell information. Actually, this kind of approach was bashfully commenced in [T^J and 
seriously pursued in [21] for what concerns the sine-Gordon theory. 

We shall distinguish two cases: 1) the repulsive regime < 77 < y (from Section 
5 this means indeed Att < b 2 < 8tc in the sine-Gordon limit); 2) the attractive case 
~ < rj < K (which covers the remainder < b 2 < An). In the first case we only have 
soliton and antisoliton excitations; in the second case we also have soliton-antisoliton 
bound states, the so-called breathers: the difference with sine-Gordon is that every- 
thing become deformed or elliptic. In general, the key ingredient shall be the primitive 
function of the integral kernel ()2.49|) . suitably normalised in the rapidity variable 9 as 

fk 9 + °° sinh 2n (rc~ 2r ?) 7r P - 2i W$ 
X (6) = / dx27rG(x) = % V 9 (K , (6.1) 

Jo w f^ 2 Sinh MK-^tt CQsh gngr „ ^ I 

which holds within the domain |Im#| < 2K. In the previous equality, we preferred to 
use the renormalised rapidity 9, which is connected to the bare rapidity x used in (12.37(1 
by the relation 

9=— x. (6.2) 
V 

In fact, the excited state version of the NLIE ()2.37() can be now written in a more 
compact form in terms of 9 by following the pattern of ^2] : 

N h N c N w N sc 

Z N {9) = NF(9) + J2x(0-h k )-J2x(0-c k )-J2xn(0-w k )-Y,Xii(0-s k ) + 

k=l k=l k=l k=l 

+ 2 / 2 ; df]G(9 - f])lm\n [l + e iZN ^ +i0) ] , (6.3) 
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for an excitation with Nh (real) holes, N c close (complex) pairs, N w wide (complex) 
pairs and N sc (complex) self-conjugated roots. 3 The reason for the appearance of the 
function xn, the second determination of x, is an important technical detail which will 
be discussed in what follows. We have also introduced a different normalisation for the 
kernel: G{9) = ^G (k^)" ^ or characterising the state, additional algebraic equations 

on Z N {6) would also be necessary [IS]; but these do not affect directly our scattering 
treatment and then are left out in this context. 



6.1 Repulsive regime: < 77 < y 

In this regime of the XYZ model, the two particle asymptotic states are spanned by 
four independent vectors, which describe respectively the soliton-soliton, the antisoliton- 
soliton, the soliton-antisoliton and the antisoliton-antisoliton excitations. Although the 
two-particle S-matrix is, as said before, the building block, we may consider, in general 
and without extra difficulties, the multi-particle states as those forming a representation 
space of the Zamolodchikov-Faddeev algebra. In fact, this is a non-commutative algebra 
with generators A{9) and A{9) satisfying 

A{9 X )A{9 2 ) = S(6i - 9 2 )A{e 2 )A(9 1 ) + S a (6 1 - 9 2 )A(9 2 )A{9 1 ) (6.4) 
A{9 1 )A{9 2 ) = S t (9 1 -9 2 )A(9 2 )A(9 1 )+S r (9 1 -9 2 )A(9 2 )A(9 l ) (6.5) 

A(e x )A{e 2 ) = s t {e 1 -9 2 )A(e 2 )A(9 1 ) + s r (9 1 -e 2 )A(9 2 )A(e 1 ) (6.6) 

A{9 X )A{9 2 ) = S(9 1 - 9 2 )A(9 2 )A(9 1 ) + S a (9 1 - 9 2 )A(9 2 )A(9 1 ) . (6.7) 

Physically, the generator A{9) (A(9)) describes a soliton (antisoliton) excitation on the 
vacuum with rapidity 9. Of course, this heuristic consideration leads immediately to 
the mentioned particle representations of this algebra: asymptotic in (out) states are 
created on the vacuum by products in which the generators are arranged with decreasing 
(increasing) rapidities. As a particular case, the scattering in the two-particle sector 
may be algebraically described by the four state basis mentioned at the beginning of 
this Subsection and the functions S(9i 2 ), S a (9i 2 ), S r (9i 2 ), S t (9i 2 ), where 9i 2 = 9\ — 9 2 , 
selected as the only non- vanishing elements of the two-particle S-matrix: 

out < A(9 l )A(9 2 )\A(9 l )A(9 2 ) > m = out < A(9 1 )A(9 2 )\A(9 1 )A(9 2 ) > m = S(9 12 ) , (6.8) 

out < A{9 1 )A{9 2 )\A{9 1 )A{9 2 ) > m = out < A{9 l )A{9 2 )\A{9 l )A{9 2 ) > in = S t (9 12 ) , (6.9) 

out < A{9 l )A{9 2 )\A{9 l )A{9 2 ) > m = out < A(9 1 )A(9 2 )\A(9 l )A(9 2 ) > m = S r (9 l2 ) , (6.10) 

out < A{9 l )A{9 2 )\A{9 l )A%) > in =out< A{9 1 )A{9 2 )\A{9 1 )A{9 2 ) > m = S a (9 l2 ) . (6.11) 



3 Generally, the roots come in complex-conjugated pairs because of the real analyticity of the count- 
ing function Z(x). A pair of complex conjugate roots is said to be close if the absolute value of their 
imaginary parts is < 2K, - i.e. the domain of the function %(6>) Itj.lfl -, wide if it is > 2K. Nevertheless, 
one root may be single if and only if it is either real or self-conjugated. In other words, they behave 
exactly like sine-Gordon roots |15j . although the fixed imaginary part of one self-conjugated root will 
be -K 2 /V 
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From these considerations, it follows that each of the following four states 



-^=(\A(9 1 )A(9 2 ) > in ±\A(9 1 )A(9 2 ) > m ) , ^=(\A(e 1 )A(9 2 ) 



±\A(9 1 )A(9 2 ) > m ) 

(6.12) 

is preserved by the scattering, i.e. it is an eigenvector of the two particle S-matrix. The 
respective eigenvalues are the four amplitudes of the scattering processes 

1 



■ out < A(9 1 )A(9 2 ) ± A(9 l )A(9 2 )\A(9 1 )A(9 2 ) ± A(9 1 )A(9 2 ) > m = S t (9 12 ) ± S r (9 12 ) 



(6.13) 



'out 



< A{9 1 )A{9 2 ) ± A{9 1 )A{9 2 )\A{9 1 )A{9 2 ) ± A{9 l )A{9 2 ) > in = S{9 12 ) ± S a ( 
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As a matter of fact, we have written down the eigenvalues of the two-particle S-matrix 
because we will read off exactly them from the IR limit (N — > oo) of the excited state 
NLIE ffOjl . 

Let us start by considering the first eigenvector 



V2 



(\A(9 1 )A(9 2 ) > in +\A(9 1 )A(e 2 ) > t 



(6.14) 



which describes the symmetric combination of a soliton and an antisoliton, and by 
computing its eigenvalue 

S+(8) = S t (8) + S r (6) , (6.15) 

with the shorter definition 6 — Q\ 2 . As Bethe Ansatz state it is built up by exciting two 
holes and a (close) complex pair of roots upon the (Fermi-Dirac) real root sea. Now, 
this configuration entails in a standard manner an excited state nonlinear integral 
equation of the form ()6.3|) with a contribution of four x{@ ~ $«) = h\, h 2 , c, c denotes 
the two holes or the two complex roots respectively). However, as far as the scattering 
factor is concerned, these four become three once 9 equals any of the hole rapidities. 
Decisive simplifications occur as N grows: the imaginary parts of the complex roots 
approach ±K respectively, since the expression (J6.1J) has exactly poles at (Im^l = 2K, 
and the real part tends to the hole middle point (as anticipated the convolution integral 
becomes negligible). Thanks to all these facts, we can identify with 9 the difference 
between the hole positions and write down the scattering amplitude as the exponential 
of the new contribution 



S + (9) = exp 



-IX 



9 --iK 
2 



9 

- + iK 



exp[ix(6>)] ■ 



(6.16) 



Upon exploiting ()6.1|) . this scattering amplitude may be manipulated into 
S+(0) = exp - } _ , . X,w-.„r (,X P 



E 

n=l 



nsinh 



2n(K-r?)?r 



|^sinh 2 "( K K; 2 ^ sin^g 



n=1 nsinh 



2n(K-r?)7r 



cosh 



(6.17) 
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or expressed in a more compact manner thanks to the function 9\\ (|2.4|) with nome 
exp (-2^): 



S + (9) 



-eKK 



T)6 



2KK' 



K' ' K' / 

exp 



011 \ 2KK' ^ l K" * h K' / 



^ zsinh 2w(K K/ 2??)7r sinfg 
^nsinh ^^ cosh 2 ^ 



K' 



(6.18) 



Let us pass on to the antisymmetric state in the same soliton-antisoliton sector 

1 



V2 

and find out its eigenvalue 



(\A(9 1 )A(9 2 ) > in -\A(9 1 )A(9 2 ) 



= S t {9) - S r {9) . 



(6.19) 



(6.20) 



It may be analogously described by a configuration with two holes and a wide pair 
of complex roots and then implies a four x contribution to the r.h.s. of the excited 
equation ()6.3j) . While N is growing the two roots tend to be a single self-conjugate 
wide root with imaginary part — K. 2 /i] and real part again fixed by the hole middle 
point. In fact, this root lies outside the strip (Im^l < 2K (here < 77 < K/2) and 
therefore we need to extend x{@) (|6.1j) into the so-called second determination, which 
is given (up to a constant) by 



Xii0) = X0)+x( 
Therefore, ()6.1|) yields explicitly 



sgn(Im#)2zK) 



Xi 1 (0) = - log 



#01 ( job 



sgn(W)^ + ^;2^ 



: K-27? .q-K->7 
' K' ' Ll K' 



exp 



-2tt 



K — 2t] 



through the entire domain of second determination 

, 2K 2 

2K < \lm9\ < 2K, 

V 



(6.21) 

sgnlm^j 
(6.22) 

(6.23) 



which is the range of variation of the imaginary part too: 2K < K. 2 /r] < 2~K 2 /r] — 2K 
(as < 7] < K/2). Therefore, with the same definition of 9 the exponential of the new 
contribution reads 



S_(0) = exp 
or equivalently if we make use of ()6.22|) 

i^yO 1701 ^ 2KK' K 7 ' K' / 







-iXn (j 




V ) . 



exp[ix(9)] , 



(6.24) 



S-(0) 



-eKK' 



^01 ("2KK 7 + i K 7 ' 2i ^ 2 ) 



exp 



isinh 2 ""^, 2 "'' sin 



K' 



KK' 



? nsinh M^L cosh^j 



. (6.25) 
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Now, we want to highlight that we have chosen the undetermined constant factor in 
the definition of x{9) and XniP) m sucn a wa Y that S + (9 = 0) = 1, S-(9 = 0) = — 1: 
this is indeed what happens for the limiting amplitudes in sine-Gordon. 

At this point, we are left with the calculation of the remaining two eigenvalues, 
S\ (9) = S(9) + S a (9), S 2 (9) = S(9) — S a (9), corresponding respectively to the two 
second eigenvectors of (j6.12j) : 



1 

7! 



(\A(9i)A(9 2 ) > m ±\A(9i)A(9 2 ) > m ) 



(6.26) 



But now the similarity with the sine-Gordon configurations clearly fails. In fact, there 
[To] the antisymmetric combination is simply realised by a Bethe state with two holes 
and hence would breed a two x{6 ~ @j) (j = 1)2) contribution in the r.h.s. of (|6.3jl . 
And in sine- Gordon theory it gives the eigenvalue of the symmetric combination as 
well. In other words, both eigenvalues coincide with the scattering amplitude soliton- 
soliton into soliton-soliton there, being the soliton-soliton into antisoliton-antisoliton 
event prevented by topological charge conservation (S a = 0). Yet, in the XYZ chain 
the situation is richer and this degeneracy removed: in physical language the channel 
soliton-soliton into antisoliton-antisoliton is here allowed and the topological charge not 
conserved but modulo 4. As a consequence, we need to go along a different route which 
shall be connected with the existence of a novelty in this scenario: the real periodicity. 
And indeed, in the elliptic case, we may also consider configurations with two holes and 
a (close or wide) complex pair in which the position of the real part of the complex pair 
undergoes a shift by ±|K' (in the x coordinate) from the middle point of the holes. 
In terms of 9 = —x, the scattering factors deriving from these new configurations are 
given by (|6.18|) and ()6.25|) (respectively for symmetric and antisymmetric state) with 9 
replaced by 9±^- (the hole contribution exp[ix(#)], which should not change, is indeed 
invariant under those shifts of 9). These simplify further after using #11(2: ± 1/2; r) = 
T9iq(z; t) and 9qi(z ±1/2; r) = #00 (z; r) respectively: 
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2»V) \^ 1 sinh g^gg sin ^ 

0oo(-^7 + ^;2z^) P ^nsinh^^cosh^ 



(6.27) 



Up to now, we found the scattering amplitudes in the basis vectors describing two 
particle states (of solitons and antisolitons) . However, in order to find a link with the 
famous Baxter's eight- vertex R- matrix (j2Hl and [2Sj), we need to describe the same 
scattering process in an equivalent way by using a different basis. In this manner, we 
will also deduce the Zamolodchikov's S-matrix of jHJ. In fact, we may introduce the real 
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doublet of particles, A\ and A 2} as 

MB) = ^[A 1 {e)+iA 2 {9)] , A(0) = -±=^(9) - iA 2 (9)) . (6.28) 

In terms of these new generators the Zamolodchikov-Faddeev algebra (J6.7)) looks as 
follows 

A 1 (9 1 )A 1 (9 2 ) = a(9 1 -9 2 )A 1 (9 2 )A 1 (9 1 ) + a a (9 1 -9 2 )A 2 (9 2 )A 2 (9 1 ) (6.29) 

A 1 (6 1 )A 2 (6 2 ) = a t (9 1 -9 2 )A 2 (9 2 )A 1 (9 1 )+a r (9 1 -9 2 )A 1 (9 2 )A 2 (9 1 ) (6.30) 

A 2 (6 1 )A 2 (6 2 ) = a{9 l -9 2 )A 2 {9 2 )A 2 {9 1 ) + a a {9 l -9 2 )A 1 {9 2 )A 1 {9 l ) (6.31) 

A 2 {9 l )A 1 {9 2 ) = a t {9 1 -9 2 )A 1 {9 2 )A 2 {9 1 )+a r {9 l -9 2 )A 2 {9 2 )A 1 {9 1 ) 1 (6.32) 

where the cx-amplitudes are related to the S-amplitudes by means of the relations 

2a = S t + S r + S + S a , 
2o" a = St + S r — S — S a , 
2a t — S + S t — S a — S r , 

2a r = S + S r - S a - S t . (6.33) 
Namely, the eigenvalues are given by 

S + = S t + S r = o + a a , 

S- = S t — S r = O t — O t , 

S\ = S + S a = o - a a , 

S 2 = S - S a = cr t + a r ■ (6.34) 

From the expressions for S + , S-, Si, S 2 (j6.18l 16.251 I6.27|) and after some lenghty 
calculation, we obtain the cx-amplitudes, 



- 1 ; P)(x- X ; Q 4 )(Q 2 x; Q^^P; P 2 )(aP; P 2 )(Q 2 ; P 2 )(P 2 Q- 2 ; P 2 ) 



(x; P)(x; Q 4 )(g 2 x-!; Q^x^Q 2 ; P 2 )(P 2 Q- 2 x; P 2 )(P; P 2 ) 2 
where, with a new definition of x which holds hereafter 4 , 

X = 6XP (^K^)' ^ = exp (~^")' P = exp ( 4 ~EC~/ ' < ^ 6 ' 36 ' ) 
and also the remaining ones 

^)_ sn(|f;4^) ffr(g )_ m (g*-ig*;4^) 



<r(g) V KK' K' K' y V KK ' ' K 



4 Because of the heavy technicality, we preferred to keep the symbol mainly used in the current 
literature at the cost of this notation abuse. 
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Srep(@) 



(6.38) 



Of course, the Jacobian elliptic function sn has nome P now and also K and k are 
respectively the first kind complete elliptic integral and the modulus corresponding to 
the same nome P. 

We are now ready to identify our new a-amplitudes with the entries of the eight- 
vertex R-matrix. Moreover, Baxter's R-matrix was used more recently in order to 
define the elliptic algebra A qe;Pe (sl(2) c ) [27J. More precisely, either the Baxter's R- 
matrix RB(x e ;q e ,p e ) (defined by (22) in |2Zj) and the R-matrix with scaled nome 
R* B (x e ; q e ,Pe) = RB(x e ;q e ,pl = p e qe 2c ) are involved in the definition of this elliptic 
algebra. Therefore the relevance of this identification, which goes through the simple 
definition 

/ a r (9) a t (6) \ 
a a {6) a{6) 
a{6) a a (6) 
V cx t (6) a r {B) ) 

As usual we shall parametrise the elliptic affine parameter x e by an exponential mapping 
of the physical rapidity 

-2 ( 2i7T?7^\ 

x e = x = exp I J . (6.39) 

Furthermore, we need to relate the deformation parameters, q e = Q and p e g~ 2 = P, to 
finalise our link 

R* B {x e ; q e ,p e )\c=i = S rep (9) . (6.40) 

Regarding the previous relation it is worth saying that the authors of |2Zj highlighted the 
impossibility to find the S-matrix of the XYZ chain in the literature, but reported a con- 
jecture due to F. Smirnov according to which it should be given by —R* B (x e ; q e ,p e )\c=i- 
In this respect, we are now in the position to unveil the mapping between the algebra 
parameters q e , p e and the physical variables of the XYZ chain, 

( 2r?7r\ / 4ttK\ 

q e = exp -— — , p e = exp — — , (6.41) 



K' J 1 V K ' / 

along with the previous relation ()6.39J) concerning the rapidity. As a check, we can 
verify that the range of parameters considered in |2Tj . 

0<Pe<qt, (6-42) 

corresponds indeed to the repulsive regime. 

As a second comparison, we want to show up the way of relating our results on XYZ 
scattering factors to the Z 4 -symmetric S-matrix obtained by A.B. Zamolodchikov 
The latter was derived as an elliptic solution of the factorization (Yang-Baxter), unitar- 
ity and analyticity conditions depending on the rapidity 6 Z and the parameters 7 and 
7'. Once we identify these with the XYZ chain rapidity and parameters, respectively, 
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in this manner 

(2mri6\ ( 'U-k9 z \ ( At]tt\ ( 4vr 2 

ex HwJ= ex H~J ' ex H"K 7 J =exp v ~ 

exp ^^^tt) = exp ^-4tt^J , (6.43) 

we obtain that our a, a a , a t and a r coincide with the homonymous quantities in [5]. 

The last due comparison is with a Takebe's work [2*5] . which apparently do not con- 
tain mention to Zamolodchikov's S-matrix. This paper contains several subtleties which 
make the verification of its Tightness quite difficult (and probably did not contribute to 
its diffusion). Nevertheless, the method is remarkably interesting and describes XYZ 
chain states directly in the N — > oo limit by converting Bethe equations into a linear 
equation for the "density" of roots. Scattering factors for excited states are then read 
off from the o(l) next-to-leading contribution to the phase shift. The latter arises after 
a complete circulation of the spatial (periodic) direction. If Takebe's parameters are 
expressed in terms of ours as 

~ 2K ' TT= 2K' " r= 2K' (fM4) 

it is rather an easy matter to reformulate our scattering factors in terms of those in 
(2.3.54) of HU: 

S+(6) = (2.3.54)/l , S_(0) = (2.3.54) /3 , -Sx(9) = (2.3.54)/2 , -S 2 (6) = (2.3.54) /4 

(6.45) 

where (2.3.54) exploits in "23] the new variables xt = — and bp — 

Finally, in order to have a check on our results, it is important to analyse what 
happens to the XYZ repulsive S-matrix in the scaling sine-Gordon limit, i.e. K' — ► 
oo, which implies K — ► 7r/2 (cfr. Section 5). In this context, thanks to the useful 
trigonometric limit 

f2rj9K A .K-rj\ sinh ^ 

7T — 2)7 

we may easily obtain the sine-Gordon values 



lim sn^-4^U Sinh ^ (6 46) 

k^co^KK'' 4 ' K' J cosh^' [ ] 



S+{6) 


= S t (6) 


+ S r (6) - 


sinh vie+ ^ 

7T-2r? c i XsG (e) 

sinh v{e -^ 




S-(0) 


= S t (9) 


- S r {6) - 


cosh v{e+ ^ 
cosh 2^ 

7T — If] 


(6.47) 




= S(6)- 


f S a (9) - 


4 e iXsG(6) 




s 2 (e) 


= S(6)- 


- S a {6) - 


+ e %o(0) 
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with the overall factor 

JxsgW _ 



exp 



dk 



sinh fc I - 



k sinh fc(£-§) cosh ^ 



Arj 



sin &6 1 



(6.48) 



describing the sine-Gordon scattering soliton-soliton— >soliton-soliton jTH]. 
6.2 Attractive regime: y < 77 < K 

As anticipated, it is really worth to complete the picture by studying the scattering 
theory in the attractive regime K/2 < rj < K: here, soliton and antisoliton excitations 
do not enjoy apart existences only, but also bound states. And we fancy to call these 
states (elliptic) breathers in analogy with the sine-Gordon locution. Although such 
bound states - whose number depends indeed on the value of r\ - shall have different 
scattering factors, for simplicity's sake we restrict our attention to the lightest one. In 
fact, we may expect that all these scattering processes should be somehow described by 
the algebraic structure of our simplest case. At any rate, their account would deserve 
a separate work [2*Hj . 

Let an asymptotic state of the lightest breather be created with rapidity 9^ i — 1,2, 
by a Zamolodchikov-Faddeev generator B(6j) with exchange relation given by a two- 
particle amplitude Sb(@i — #2): 



B(9 1 )B(9 2 ) = S 



9o)B(9?)B( 



(6.49) 



In our set-up Sb{9% — 62) derives from considering the thermodynamic limit of the 
NLIE describing the excitation of two lightest (elliptic) breathers. This configuration 
with two breathers B{6i) with rapidities 0j corresponds to adding up, to a sea of real 
roots, two self-conjugate roots with real part 9i respectively, and imaginary part K. 2 /r) 
(limiting value). Then, the NLIE ()2.37|) takes on the form (discarding the convolution 
term as usual) 

Z{9) = NF{9) - xii (0 -9,- 1—) - xii (§ -9 2 - 1—) , (6.50) 



when |Im#| < 2K 2 /r/ — 2K is in the first analyticity strip. Nevertheless, the imaginary 
part of a self-conjugate root — > 2K 2 /r/ — 2K turns out to be outside the first analyt- 
icity strip and therefore we had to use the second determination of the function x(^)j 
i.e. 

Xn(0) = X (0) - X (0 - 2i(K 2 / V - K)sgnlm#) , (6.51) 

which holds as the imaginary part |Im#| > 2K 2 jr\ — 2K lies outside the first analyticity 
strip. Implementing now the definition (|6.1jl . it follows a compact formula for the 
functions Xii(@) i n Q6.50J) 

(x-'P-'Q 2 ; Q 4 )(xQ 2 ; Q^x^Q 2 ; Q 4 )(xPQ 2 ; Q 4 ) 



XlI (9) =i\n 



-1. 



Q A )(xP; Q±)(x- l P- l Q^ Q 4 )(a;Q 4 ; Q 4 ) 



(6.52) 
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where P and Q are defined in (|6.36|) and 



[2mr]d\ 



(6.53) 



As second and last step, we compute Z{6) at one of the self-conjugate roots, for instance 



e = 0! + i- 



K 5 



(6.54) 



and expect that the terms involving \ii furnish % In Sb{Q\—Q'i)- And again the imaginary 
part of a self-conjugate root (|6.54|) — > 2K 2 /r^ — 2K lies outside the first analyticity 
strip of Z(9) itself: therefore we need to consider its second determination as well. In 
other words, we must introduce the second determination of Xn{®)i Xn(@)u : 



Zn (e x + i^j = NF (ex + i^J - XiMn - Xii(0i - 02)11 ■ 

We can calculate it by reiterating ()6.51|) with now lm9 > and simply obtain 

(x^P-'Q 2 ; Q 4 ){x~ l P- Q^xP-'Q 4 ; Q 4 )(xPQ 2 ; Q 4 ) 



(6.55) 



Xi 1 (0) 
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i In 



-x- 



(xP-^Q 2 ; Q 4 )(xP; Q^x^P^Q 4 ; Q^x-iPQ 2 ; Q 4 ) 
where we have eventually identified 

9 = 61-62, x = exp f • 



(6.56) 



(6.57) 



In conclusion, the scattering factor between the lightest elliptic breathers reads in this 
XYZ notation 



Sb{6) = -exp [ixii{0)n] exp [ixn(6)ii] = 
1 (xP-'Q 2 ; Q 4 ){xP; Q^jx^p-'Q 4 ; Q 4 ){x~ l PQ 2 - Q 4 ) 
x {x-ip-iQ 2 ; g 4 )(a;- 1 P; Q^xP^Q 4 - Q 4 )(xPQ 2 ; Q 4 ) 



(6.58) 



with the above 6 and x. An overall minus sign has been permitted thanks to the 
definition of xn U P to a constant and in order to reproduce Sb(0) = —1 as in the 
sine- Gordon field theory. 

Despite the cumbersome calculation and the different origin, this scattering factor 
coincides exactly with the structure function of the Deformed Virasoro Algebra Vir Pv ^ qv 
introduced by Shiraishi et al. |4J . The Deformed Virasoro Algebra (DVA) is an asso- 
ciative algebra generated by the modes T n of the current T(z) = ^2 n T n z~ n , satisfying 
the relation 



f(w/z)T{z)T(w) - f(z/w)T{w)T{z) 



(1 - -p v q v 1 ) 
1-Pv 



8 



p v w 



5[ — 

\PvZ / 



(6.59) 
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where p v and q v are complex parameters and 

(i-^)(i 



f(x v )= exp 



n=l 



1 +p" n 



(6.60) 



Equality ()6.59|) is to be interpreted as an equality between formal power series, but, 
as shown in Appendix B.l, f(x v ) can be analytically continued to the whole complex 
plane and this allows us to recast relation (|6.59jl in the braiding form 

T(z)T(w) = Y(z/w)T(w)T(z) , (6.61) 

where 12 211221 2 

y/ \ 1 (x v q v p Vl p v ) (x v q v , p v ) (x v q v p v ,p v )(x v QvPv,p v ) ,„ 

x v \x v q v Pv,P v ){% v Qv)Pv)\ x vQ v P v ] P v )\ x vQvPv] P v ) 

is the structure function of the DVA. While writing this formula, we have borne in mind 

the convenient infinite product notation 

00 

(x-a) = JJ(1 -xa s ) , (6.63) 

introduced in (|2.5j) . Eventually, if we suppose the mapping between spin chain variables 
and algebra parameters 

/ 2mri6\ ~ ( 4ttK\ „ / Ann 

x v = x = exp f — iciC 7 " J ' Qv = = 6XP V — K 7 " J ' Pv = = GXP V "I 77 

(6.64) 

we can see immediately that 

S B (6)=Y(x v ). (6.65) 

Therefore, we have proved that the Deformed Virasoro Algebra defines a current, T(z), 
which closes the Zamolodchikov-Faddeev algebra for the fundamental scalar excitation 
of the XYZ model. 

It is possible to rewrite Sb{8) in a different form, which is useful if we want to 
compare this factor to analogous ones proposed in the literature. We follow the idea of 
Lukyanov jB], write (j6.65l6.62j) in terms of theta-functions, perform a modular trans- 
formation on the elliptic nome and eventually obtain (details in Appendix B.2) 

ft ( *3 I K . iK.' \ n ( J0 K . iK' 

S B {8) = 4K ^ * V l / K 2n,AT '( . (6.66) 

a I i6 K . iK' \ a ( iO K.iKM 

° n \4:K 27?' 4r? J Ul ° \4K ^ 2ij ' 4r t J 

This is a useful form as it allows us to perform easily the limit K' — > 00 which describes 
- after shifting (Section 5) - the sine-Gordon model. Indeed, we obtain the first breather 
scattering factor [TH] 

sinh 9 + i sin tt£ 

hm S B {0) = — — — , (6.67) 

K'-+oo sinh 8 — 1 sin 7r£ 
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with the usual parameter 

«-5j-l. (6.68) 

Moreover, it follows from (|5.13|) that £ carries on the dependence on the coupling 
constant b of the sine-Gordon Lagrangian: 

It is an easy job, now, to relate ()6.66|) to the scattering factor Smp{P) proposed by 
Mussardo and Penati in jjj. We start from the form (6) in [7j for Smp(P) and after 
some manipulations (details in Appendix B.3), we can conclude that this expression 
(6) equals our ()6.66|) . namely Smp(P) = Sb(0), provided these identifications are taken 
into account (parameters of [7j are on the l.h.s.) 

ttK' _ K n6 

Tmp — — — , (Imp — 1 , Pmp — ~ ttf? ■ (6.70) 

2rj rj 2K 

As an obvious but intriguing byproduct we want to underline the coincidence of the 
Mussardo-Penati's scattering factor and the structure function of the DVA. 

Remark: The parameters x v , q v and p v of the attractive regime (and of the DVA 
structure function) inherit a simple expression in terms of the parameters x e , q e and p e 
of the repulsive regime (and of the R- matrix of the elliptic algebra A q ^ Pe (sl(2) c )): 

X v =X 2 e , Pv = ql, Qv=Pe- (6.71) 

Although in a different landscape, this seems to be exactly the connection between 
Vir Pv>qv and A qe , Pe (sl{2) c ) proved in 



7 Some conclusions, many prospects 

Starting from the (elliptic) Bethe Ansatz equations of the spin 1/2-XYZ chain on a 
circumference (cfr. for instance |2*lj). we have written a Non-Linear Integral Equation 
describing describing a generic state (either the vacuum or an excited state) of the 
model in the disordered regime. Maintaining the size finite, we have studied two dif- 
ferent limits: the usual trigonometric limit which gives the spin 1/2- XXZ chain and a 
double scaling limit which turns out to describe the sine-Gordon field theory on a cylin- 
der. The latter furnishes the generalization of the infinite length scaling limit of 2j and 
[3j in case of finite volume. Moreover, it has suggested us the heuristic idea that send- 
ing to infinity the size of the 1/2-XYZ would infer an elliptic scattering theory, which 
naturally inherits an elliptic deformation of all the sine-Gordon structures. In fact, any 

5 As a conjecture this liaison is already in PU] and somehow in . 
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elementary excitation on the vacuum in the repulsive regime can be completely charac- 
terised, through a NLIE, by new terms which give rise to the corresponding scattering 
amplitude. Therefore, the elliptic deformation of the soliton/antisoliton sine-Gordon 
S-matrix derives from the finite size procedure and the re-expression of this matrix as 
the Baxter elliptic R-matrix proves what is called in [2Zj the Smirnov's conjecture, a 
remarkable connection between representation theory of an elliptic algebra and a scat- 
tering S-matrix. Very likely the conjecture may also be extended to more complicated 
elliptic algebras, being its formulation purely algebraic. 

Moreover, we have studied the lightest bound state, i.e. the elliptic deformation 
of the first sine-Gordon breather, and found that its scattering factor coincides with a 
proposal by Lukyanov jHj and, in a different context, by Mussardo-Penati |7j 6 , provided 
a mapping between each set of coupling constants and scattering rapidities is given. 
Actually, the most important result of this calculation is that the scattering factor 
coming from the spin chain shows manifestly its identity with the braiding factor of the 
Deformed Virasoro Algebra by Shiraishi-Kubo-Awata-Odake jl] and let us suppose that 
each heavier breather may produce a braiding algebra with a similar structure (possibly 
like in [HI], but without one of the squares in the structure function). In this respect, we 
would like to clarify this point in a ongoing publication |2H], since we reckon the basic 
DVA of |P as responsible somehow for the algebraic structure and the mass quantisation 
of the other breathers as well. This conjecture relies upon the identification of the first 
breather as the fundamental scalar particle and upon the vertex operator construction 
of DVA given in [22] ■ We would like to conclude this part about the relevance of the 
DVA by highlighting how its field theory limit in sine-Gordon should be the quantum 
version of the non-local symmetry geometrically constructed in [33] . But the field theory 
inheritance really needs an ad hoc analysis in a separate paper, though the scenario is 
becoming clearer thanks to these spin chain developments. Nevertheless, it is also worth 
analysing the whole information obtainable on the lattice system through at least two 
routes: on the one hand, pretty much is now known about the representation theory 
of the Deformed Virasoro Algebra (cfr. j3U] as a review work); on the other hand, the 
form factors postulates authorised Mussardo and Penati to work out an entire series of 
form factors of not better identified "fields" [7} (where also some speculations on the 
nature of the fields is brought forward in the light of conformal fields structure). By 
now, a more precise reading of these form factors becomes possible and a field theory 
interpretation plausible. 

In conclusion, a field theory description of the spin chain in the thermodynamic 
limit, also conveying and re-interpreting previous results, shall be pursued by means 
of fermionisation and renormalisation group techniques and would contribute to the 
clarity of the landscape. Towards this aim, valuable results from J2U] might be efficiently 
extracted. 

6 The coincidence of these two was apparently unnoticed beforehand. 
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Appendix A 



In this appendix we want to compute the Fourier coefficient of the function 

d 



dx 



(x + ia,rj) , a G 



(Al) 



From the definition (j2.17|) of </> it follows that we have to consider the expression 
—(j)(x + ia,ri) = i-^j- In H(r) + a — ix) — i— In H(r) — a + ix) + 

VAjJU KajiAj U/tA/ 

+ i— In 0(77 + a — ix) — i— In Q(r] — a + ix) . (A. 2) 

dx dx 



Using formula (8.199.1) of [22] we have: 
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Since — ^ < a; < we can express the denominators as power series: 
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Now, we are ready to perform the integration 
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About the terms depending on the cotangents we have that 
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On the other hand, the remaining terms yield 
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Summing these two expressions we get the result for (jA.4|) . However, before doing such 
a sum, it is convenient to compute the contribution coming from the function: 
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Using formula 8.199.4 of [22! we § e t 
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Since -y < i < y we can express as power series the denominators: 
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With the help of ()A.9|) we can perform the integration involved in (|A.7I 
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Now, summing (jA.6|) and (jA.lOjl . we get 
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We notice that this expression exactly cancels the term in (|A.5|) proportional to (— l) n . 
Therefore, we get the following result for the Fourier coefficient of -^(j)(x + ia, f|), a 6 I: 
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The sum over S can be performed using formulae 1.445.1,2 of |22j- The final result 
depends on the range of values of a and rj. We remark that we can remove the restriction 
to a real, since the imaginary part of a can be easily implemented in the final formula, 
its effect being a phase. Because of the periodicity property 0'(x + 2iK) = 0'(x), we 
can restrict a to the interval — rj < Rea < 2K — rj. We remember also that < rj < K 
(I2~TTJ) . We get that, for a G C: 
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r] < Rea < 2K - 77 . 



On the other hand, formula (jA.12j) holds also if we replace r\ with 27/. Since now 
< 2?7 < 2K, in order to compute the sums over S we distinguish the following cases: 
• < 7] < K/2: 
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277 < R e « < 2K - 2t/ ; 



K/2 < 7/ < K: 
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Appendix B: DVA as scattering algebra: a unifying 
outlook 



B.l XYZ and DVA: an equivalence 

In the definition of the Deformed Virasoro Algebra (DVA), the main ingredient is the 
formal expression 



f(x v )= exp 



n=l 



n 



:b.v 



The series contained in (jB.ljl is convergent when \x v \ < 1 and if, for instance, the 
parameters q v , p v satisfy \p v \ < 1, \q v \ < 1 and Ip^g^ 1 ] < 1. In such a case we can sum 
up the series and define a meromorphic function in the complex plane as its analytic 
continuation. Now, we want to determine such a function. 

We observe that, since \p v \ < 1, it is possible to write the denominator in (jB.lj) as 
a geometric series: 



f{x v ) = exp 



/ i 
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s=0 



(B.2) 



36 



Since \q v \ < 1 and *| < 1, we can exchange the series in (|B.2|) and perform the 
summation over n in a suitable range of values for x v by using the formula 



00 z n 1 

J2 — = log - — 

z — ' n 1 — 

n=l 



l^l < 1. 



(B-3) 



So we obtain 

f(x v ) = exp 



00 / 1 1 1 1 
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(B.4) 

This expression can be rearranged by separating the odd and even s contributions and 
by some simplifications 
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Using notation ()2.5|) for infinite products, namely 
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we can write down the compact expression 
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1 x v (x v q v p^ p^) (x v q v p v , p^) 
Now, we define a meromorphic function, the so-called structure function of the DVA, 



as 
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fix' 1 ) x v (x v 1 q v 1 p v ;pl)(x v 1 q v ;pl)(x v q v 1 pl;pl)(x v q v p v ;pl) 
This function enters the exchange algebra 

T(z)T(w) = Y(z/w)T(w)T(z) , 

which is of course equivalent to the DVA. 



(B.8) 



;b.9) 



B.2 Breather factor Sb(0) in Lukyanov's form 

We want to elaborate the DVA expression (|6.651 I6.62|) for Sb(0) following an idea by 
Lukyanov (6,). First, keeping in mind (|6.64j) . we collect the infinite products of ()6.62j) 
within theta-functions of nome p v = exp 



a ( V e I o • K , 4i?7\ n ( t}0 o • K . 
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Then, we use the modular transformation 

On (pr,i%) _; 0u Hf^f) mm 

n ( <K'\ a I K -K\ ' l- D - ii J 

01 {a;i ir ) Oio{-i w a;i w ) 

which connects theta-functions with nome p = exp (— vr^) (on the l.h.s.) with theta- 
functions with nome p' = exp (— ^j^r) (on the r.h.s.). This allows us to re-express Sb(9) 
as follows 

S B (9) = 4K ^ 4?? ° r (B.12) 

which is (Ifj.fjfjjl of the main text. Upon comparing this relation with (14) of [Hj, we must 
identify Lukyanov's variables (on the left) and ours in this manner: 

= ^'j- 1 - '^(-w)- (B13) 



B.3 Breather factor Sb(Q) in Mussardo-Penati's form 

Here we want to show that Mussardo-Penati scattering factor Smp{/3) actually coincides 
with Lukyanov's one and therefore with our Sb{9). We start from expression (6) of [Zj 
for Smp{(3) and simply re-formulate it by expressing the Jacobian elliptic functions in 
terms of theta-functions: 

C (a\ _ 11 I 2in ' 7T I 10 I 2in ' tt [ 01 I 2in ' tt ) 00 I 2m ' tt J /r> 1 A\ 

As stressed the elliptic nome of these theta-functions is e~ T . Now, we use the simple 
identities (deriving from the definitions), 
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9 11 (u;t)9 01 (u;t) = 9 11 (u;r/2) K ' ' e~r , (B.15) 

\e , e j 

6 10 (u;t)6 00 (u;t) = 6 10 (u;r/2) \ ; , (B.16) 
in order to express Smp(P) in the final form 

C ( R\ — 11 V 2ct ' 2tt) Ul ° \ 2iw ' 2ti-/ m 171 

y H k~2Sr' 2^J y 10 2^) 

Indeed, we notice that this coincides with Sb{6) . provided we link Mussardo- 

Penati's variables (on the l.h.s.) with ours (on the r.h.s.) according to 

ttK' K it6 

Tmp — — — , clmp — 1 , pmp — — tttt • (B.18) 
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